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Abstract. - A relationship between the measurable monomer-monomer structure factor, 
and the centre-of-mass (CM) structure factor of dilute or semi-dilute polymer solutions is 
derived from Ornstein-Zernike relations within the "polymer reference interaction site model" 
(PRISM) formalism, by considering the CM of each polymer as an auxiliary site and neglecting 
direct correlations between the latter and the CM and monomers of neighbouring polymers. 
The predictions agree well with Monte Carlo data for self-avoiding walk polymers, and are 
considerably more accurate than the predictions of simple factorization approximations. 



There have recently been a number of attempts to represent polymer coils in solution as 
systems of soft, penetrable particles of fixed or variable ^ shape, and size proportional 
to the radius of gyration Rg, which interact via effective pair forces obtained by averaging 
over individual monomer degrees of freedom. A similar approach has proved very successful 
in the description of star polymer solutions [Q . The advantage of this reductionist strategy 
is that the effective pair forces act only between single interaction sites within each polymer, 
rather than between the large number of individual monomers or Kuhn segments belonging to 
the interpenetrating coils. A natural (but by no means unique) choice of the single interaction 
site is the centre-of-mass (CM) of the polymers. This coarse-graining leads to an enormous 
reduction (by a factor equal to the number v of monomers or segments in each polymer) in 
the number of interacting degrees of freedom, thus allowing, inter alia, efficient simulation of 
large scale phenomena involving many polymers. 

The price to pay is that the effective interactions are state-dependent, i.e. the pair po- 
tentials resulting from the coarse-graining procedure depend in general on polymer concen- 
tration, temperature, and the degree of polymerization, although the latter dependence is 
negligible for sufficiently large v such that the scaling regime is reached in practice. In ref. 
the state-dependent effective pair potential between the CM of non-intersecting self-avoiding 
walk (SAW) polymers in dilute and semi-dilute solutions was derived from a very accurate 
Statistical Mechanics inversion procedure of the CM pair distribution function gcci^). The lat- 
ter had to be determined for each polymer concentration from fully microscopic Monte Carlo 
(MC) simulations of samples involving hundreds of SAW polymers. This "exact" procedure is 
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obviously very computer intensive, and in a certain sense defeats the original purpose, which is 
to replace a detailed, monomer level description of polymer coils by a projected representation 
involving only their CM's. 

To overcome this methodological bottleneck, one is naturally led to the use of approxi- 
mate theories of intermolecular correlations between polymers. One obvious candidate is the 
"polymer reference interaction site model" (PRISM), which provides a successful theory of 
such correlations, particularly so in the melt, where polymers are known to behave as Gaus- 
sian coils jsj. However PRISM provides only monomer- monomer pair distribution functions 
gramif), and an accurate procedure is needed to extract the CM-CM pair distribution func- 
tion gcc{i") from a knowledge of gmmi'i')- In this Letter the PRISM formalism is extended to 
derive such a relationship, which also involves the form factors (or internal structure factors) 
of individual polymers. The relation is tested against MC simulation data, and turns out to 
be vastly superior to previous proposals 

Consider a monodisperse solution of N linear polymers with v monomers or segments each, 
in a volume V\ the polymer number density \s p — N/V. Let R* (1 < i < N) be the CM 
position vector of the ith polymer, and = R' + uj^ {I < a < the positions of the ly 
monomers in that chain. The Fourier components of the monomer and CM density operators 
are: 

N V 
i—l a— 1 

i=l 



It proves also convenient to define the Fourier components of the intramolecular monomer 
density of the iih polymer: 

u 
a=l 

in terms of which one may express the following form factors, which characterize the internal 
structure of each individual polymer coil, in the presence of all surrounding polymers: 

w™™(g) = < PqP'Lq > , (2a) 
Wcm(g) =< Pq > • (2b) 

The brackets denote canonical averages weighted by the Boltzmann factor involving the total 
interaction energy between all monomers on all polymers, and accounting for the connectivity 
constraints. The resulting form factors are independent of the index i, because all polymers 
are equivalent, and depend only on the modulus q — |q| of the wave vector due to rotational 
invariance (the polymer solutions are isotropic). 

The total monomer structure factor, as measured e.g. by coherent neutron or light scat- 
tering experiments, is defined by 

I j a p 
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and naturally splits into intramolecular (i = j) and intermolecular (i ^ j) contributions: 

Smmiq) — S^m^il) + S'^m'^ il) ~ ^ + 'S'mnT('J') • (4) 

The intermolecular contribution may be rewritten as: 

«'(<z) - ^ E E < ^-"^-"^'f^f^. > ■ (5) 

i j^i 

A common decoupling approximation is to assume that the intramolecular conformations of 
any two polymers are independent of each other and of the mutual positions of their CM, so 
that the statistical average in ^ factorizes according to: 



S'l^'-id) - ^ E E < X Pli >< >= ^^^[Scc{q) - 1] , (6) 



^C7n{q) 

^ e - • P'-a ^= 

i j^i 

where Scc{q) is the CM structure factor: 

Scciq) = ^ < >= 1 + phcM , (7) 

and hcc{q) is the Fourier transform of the CM-CM pair correlation function hcdr) = gcc{r) — 1- 
Substituting eqs. (||) and @) into eq. (0), one arrives at 

Smm[q) = h phrnmiq) = h p 5 hcc[q) , (8) 

which leads to the desired approximate relation between the Fourier transforms of the monomer- 
monomer and CM-CM correlation functions: 

hcc{q) = r^^mm(g) ■ (9) 

Equation is the approximation proposed by Koyama ||^ and is formally identical to the 
so-called "free-rotation" approximation for rigid molecules; it requires a knowledge of the form 
factor LJcm{q) defined by eq. (pb|). 

An alternative approximation is based on the rigid particle assumption used by Pago- 
nabarraga and Cates |^ , according to which intramolecular conformations and CM positions 
are independent, and intramolecular conformations of different polymers are correlated as if 
they belonged to the same polymer ("rigid particle" assumption), i.e.: 



-Scciq)^ < PqPlq > 
^mm{q) 



leading to 



Scciq) , (10) 



hcciq) = f-rhmmiq) ■ (11) 

l^mmKq) 
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The relations (|^) and (^) are clearly based on uncontrolled factorization approximations; 
they will be tested below against "exact" simulation data, and compared to the relation which 
we now set out to establish within the framework of the PRISM theory |^ . PRISM is based 
on the assumption that all correlation functions between monomers are independent of their 
positions along the chain, i.e. end effects are neglected, which is true in the scaling limit 
u — » oo. This unique correlation function hmmiq) is related to a unique monomer- monomer 
direct correlation function c„im{q) by the PRISM Orstein-Zcrnike (OZ) relation: 

hmmiq) = ^mm{q)Cmm{q)[^m.mi{q) + Pvhmm{q)] ■ (12) 

The key idea now is to consider the CM of each polymer as an additional non-interacting 
site, which is linked to the monomer position vectors by the defining "connectivity" constraint 



R 

This trick of introducing auxiliary sites to compute special correlation functions dates back to 
the early days of the "reference interaction site model" (RISM) theory (from which PRISM 
is an extension) [Q and has been first proposed by Chandler Since the CM auxiliary site 
does not interact with any of the v segments, it must clearly be treated separately; thus, each 
polymer now has two species of sites, namely the CM and the v equivalent interaction sites 
associated with the physical segments. 

The single OZ relation (|lj) is now replaced by a 2 x 2 matrix of OZ relations. The latter are 
further simplified by the plausible assumption that the direct correlation functions between 
the CM of one polymer, and the CM as well as the interaction sites of the other polymers are 
identically zero, i.e. 

Ccc(g) = 0, (14a) 

Cmc{q) = CcmXq) = • (14b) 

The four coupled OZ relations are then given by (|2|), which is not modified by the presence 
of the auxiliary, non-interacting site , together with 

hcm{q) ^ (^cmiq)Cmm{q)[uJmmiq) + Pl'hmm{q)] , (15a) 
hmc{q) = LUmm{q)Cmm{q)[uJmc{q) + pyhmc{q)\ , (15b) 
hcc{q) = U}cm[q)Cmm{q)[^^cm{q) + pvhcm{q)] ■ (15c) 

These equations immediately lead to the desired relation: 



{qf 



hcc{q) = ,^\2 ^rnm{q) , (16) 



which is independent of any specific closure relation, except for the assumptions (p^. 

Equation ( p^ ) provides the desired link between the monomer-monomer correlation func- 
tion, which may be calculated from PRISM theory, and the CM-CM correlation function which 
is the basic input into the coarse-graining scheme, since the state dependent effective pair po- 
tential between the centres of mass follows directly from standard inversion procedures |l|. 
To extract hcc{q) from hmmiq) ^^^o requires a knowledge of the form factors uJmmiq) and 
ij^cmiq)- These are known explicitly for Gaussian chains ||l^, and PRISM theory generally 
uses the Gaussian chain form factor in the OZ relation (p^, which is strictly valid only in the 
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Fig. 1 - Density dependence of the ratio uicm{qY /o^mmisif' as a function of the wave number scaled 
by the density-dependent radius of gyration Rg{p), as obtained by MC simulations of f = 500 SAW 
polymers on a simple cubic lattice. For comparison, the corresponding function for Gaussian chains 
is also shown. 



melt. In polymer solutions, we expect the form factors to depend on polymer concentration. 
To investigate this dependence, we have carried out MC simulations of = 500 SAW poly- 
mers on a simple cubic lattice for reduced densities p/p* = 0, 0.29, 1.16, and 4.63, where the 
overlapping density p* is defined through the relation 4:Trp*Rg{p = 0)/3 = 1. Results for the 
ratio uJcmiq)^ /^mm{(lY 1 relevant for the mapping (|l6|) of hmm{q) onto hcc{q), are shown in 
fig. |l|. They are compared with the corresponding ratio for Gaussian chains, which is found to 
qualitatively reproduce the shape of the curves obtained by computer simulation. When the 
wave number q is scaled with the radius of gyration Rg (p) appropriate for each density, which 
decreases from 16.8 to 14.8 lattice spacings as the density increases, the low density results 
{p/p* ^ 1) E^re seen to fall practically on a single master curve. The results at p/p* = 4.63, 
which is well into the semi-dilute regime, deviate from this master curve, towards the result 
for a Gaussian chain, as one might expect. 

Having determined the "exact" form factors, we are now in a position to test the relations 
between hmm{q) and hcc{q), as predicted by the factorization approximations (||) and (|lT|), 
and by the PRISM-based equation (^6|). For the test to be meaningful, we have used as 
input the "exact" (rather than PRISM-generated) /imm(g). Figure § shows the MC results 
for Scc{q), together with the results obtained by combining eq. (|^) with each of the three 
approximations (^, (1^) and (|l^), using the "exact" form factors, at p/p* — 1.16. In the case 
of approximation (^), results based on the Gaussian chain form factors are also shown. The 
agreement between the "exact" structure factor and the predictions from eq. (^6|) is seen to be 
excellent at small g; some deviations are seen at intermediate wave numbers qRg ~ 4, where 
the "exact" Scc{q) has a small maximum. Similar agreement is found at the other densities, 
as shown in fig. ^ for p/p* = 4.63. The decoupling approximations and ( [Til ) are seen to 
fail. Interestingly, approximation (^6|) combined with the analytic Gaussian form factors also 
yields good agreement. 

The success of eq. (^6|) in obtaining the CM correlations from the monomer correla- 
tions suggests that the inverse route could also be profitable. This may be expected to be 
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Fig. 2 - CM structure factor Scc{q) for v — 500 SAW polymers on a simple cubic lattice at a reduced 
density p/p* = 1.16. The "exact" structure factor is compared with the various approximations 
proposed in the text using the "exact" monomer correlation function /imm(?) and form factors ujcm{q) 
and tJmm(g). For comparison, in the case of approximation (|l^), results based on the Gaussian chain 
form factors are also shown. 

Fig. 3 - Same as fig. ti, at a reduced density p/p* — 4.63. 



more difficult because one appears to move from less information (CM) to more information 
(monomers). In fact, attempts to invert eq. (^6|) result in problems of a rather general kind 
because tJcm('z), which can be zero, is now in the denominator, leading to spurious singulari- 
ties. LOcmiQ) for Gaussian coils does not cross zero, so that it may be used more profitably to 
obtain Smmio) from Scc{q), but the results are still not as good as what is seen in figs. |^ and 
^. We found empirically that the inverse of eq. (|l6|) still works well for qRg < 2, but for larger 
q deviations occur. While the other two approximations do not show the spurious singularity, 
they do show other spurious effects. 

To conclude, we have proposed a new relation between the monomer-monomer and CM- 
CM correlation functions of polymers based on an extension of PRISM theory. It has been 
tested for polymers in good solvent and it is found to be quantitatively accurate except for 
discrepancies at intermediate wave numbers where hcdq) changes sign and exhibits a small 
maximum, while hmmio) is a monotonically increasing function of q. This discrepancies must 
be traced back to assumptions (|l^, which are the only approximations involved here (apart 
from the neglect of end effects in PRISM, that is known to have minor consequences). The 
factorization approximations (^ and ( pi]), on the other hand, appear to perform less well. 
It should be added that approximation (|ll[) was derived for use in the melt, a regime which 
has not been tested in this paper, and where the connection between CM and monomer 
correlations may be different. 

The present analysis stresses the need for developing a reliable theory for the state- 
dependent form factors iOmmiq) and Wcm(9), in particular in the dilute and semi-dilute solution 
regimes that have been investigated here. In this respect, it is worth stressing that Wcm((?) 
has hitherto attracted very little attention probably because it is not of direct experi- 
mental relevance. Our MC results point to a near universality of the ratio ojcmio) / ^mmifl) as 
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a function of qRg in the dilute regime. Accurate form factors in the dilute and semi-dilute 
regime would also provide a crucial input into PRISM calculations of hmm{q)- Work along 
these lines is in progress. 

Various extensions of this work are currently under investigation. Firstly, using the two- 
component PRISM OZ relation it is easy to extend the present formalism to colloid- 
polymer mixtures and simple inhomogeneous situations (polymers near a wall for instance), 
that have recently attracted much attention. Secondly, the relation ( ]T6| ) does not rely on any 
assumption on the monomer interactions, so it would be interesting to test its relevance for 
other classes of polymers such as linear polymers in 9 or poor solvents, branched polymers, 
polyelectrolytes, etc. This could be very helpful in interpreting experiments, which typically 
only have access to monomer- monomer correlations. 

* * * 
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